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INTRODUCTION 

Anisotropic  elastic  materials  are  interesting  materials.  A  simple 
tensile  stress  applied  to  the  material  produces  not  only  an  extensional 
strain  but  also  a  shear  strain.  Likewise,  a  pure  shear  stress  applied  to 
the  material  produces  a  shear  strain  and  an  extensional  (or  compression) 
strain.  Therefore  a  loading  which  is  symmetric  (or  antisymmetric)  with 
the  xr-axis,  say,  in  general  does  not  produce  a  deformation  which  is 
symmetric  (or  antisymmetric)  with  the  xr-axis.  There  are  surprises  in 
which  anisotropic  materials  behave  like  isotropic  materials.  These  will 
be  pointed  out  in  the  paper. 

In  contrast  to  isotropic  elastic  materials  which  have  two  elastic 
constants,  anisotropic  elastic  materials  may  have  as  many  as  twenty  one 
elastic  constants.  When  two-dimen.  onal  deformations  are  considered, 
the  analysis  still  requires  fifteen  elastic  constants.  In  view  of  this,  there 
is  a  wide  spread  and  justifiable  misconception  that  the  analysis  of 
anisotropic  elastic  materials  is  much  more  complicated  than  that  of 
isotropic  elastic  materials.  This  is  not  necessarily  true  if  one  employs 
the  Stroh  formalism.  With  the  Stroh  formalism  the  solutions  to 
anisotropic  elasticity  problems  are  in  most  cases  simpler  than  those  for 
isotropic  elasticity  problems.  The  reason  is  simply  that  isotropic 
materials  are  more  than  a  special  case  of  anisotropic  materials.  They  are 
mathematically  degenerate  materials. 

Much  progress  has  been  made  since  Stroh’s  two  pioneering 
papers  appeared  in  1958  and  1962  [1,21.  We  will  point  out  in  the  paper 
the  integral  formalism  of  Barnett— Lotne  [3]  which  allows  us  to  compute 
three  Barnett-Lothe  tensors  S,  H  and  L,  which  are  real,  without  finding 
the  Stroh  eigenvalues  p  and  the  associated  eigenvectors  a,  b,  which  are 
complex.  We  will  also  point  out  some  identities  which  enable  us  to 
convert  certain  combinations  of  p,  a  and  b  to  S,  H,  L  and  other  real 
quantities.  Owing  to  these  identities,  several  existing  complex  form 
solutions  are  simplified  to  real  form  solutions  and  solutions  are  obtained 
for  some  heretofore  unsolved  problems.  As  a  result,  many  physically 
interesting  and  unexpected  phenomena,  which  have  been  shrouded  in  the 
complex  form  solutions,  have  been  discovered  recently.  Most  of  the 
unexpected  phenomena  defy  an  intuitive  explanation. 
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The  Ctroh  formalism  is  not  only  mathematically  elegant  and 
technically  powerful,  but  some  of  its  mathematical  quantities  such  as  the 
eigenvalues  p  and  the  eigenvectors  a  and  b  have  physical  meanings.  The 
mathematical  structure  of  S,  H  and  L  provides  us  a  rare  insight  into  the 
relations  between  anisotropic  and  isotropic  materials. 

For  isotopic  materials  the  in— plane  displacement  and  the 
antiplane  displacement  are  uncoupled.  The  in— plane  displacement  (ut, 
uj)  and  the  associated  surface  traction  vector  on  any  boundary  V  are 
polarized  on  the  (xi,  x3)  plane  while  the  antiplane  displacement  us  and 
the  associated  surface  traction  are  polarized  along  the  X3  axis.  For 
general  anisotropic  materials  under  the  assumption  of  two-dimensional 
deformations,  the  U3  component  is  in  general  non-zero  and  cannot  be 
uncoupled  from  the  in-plane  displacements  ui,  uj.  This  does  not  mean 
that  there  are  no  planes  or  axes  on  which  the  displacement  and  the 
surface  traction  are  polarized.  There  are,  as  we  will  show,  oblique  planes 
and  axes  on  which  the  displacement  and  the  surface  traction  are 
polarized. 

The  synergism  of  mathematics  and  mechanics  appears  to  work 
very  well  for  anisotropic  elastic  materials.  Examples  presented  in  the 
paper  illustrate  that  this  is  indeed  the  case. 


1.  THE  STROH  FORMALISM.  In  a  fixed  rectangular  coordinate 
system  X*  (i  =  l,  2,  3)  let  ui,  <r\ j  be,  respectively,  the  displacement  and 
stress  in  an  anisotropic  elastic  material.  The  stress  strain  laws  and  the 
equations  of  equilibrium  are 


=  Cijksuk,s  ' 

(1.1) 

Cijksuk,sj  =  °’ 

(1.2) 

where  a  comma  stands  for  differentiation,  repeated  indices  imply 
summation  and  Cjjks  are  the  elasticity  constants  which  are  assumed  to 
possess  the  symmetry  property 

Cijks  =  Cjiks  =  Cijsk  =  Cksij  * 

For  two-dimensional  deformations  in  which  Ui  (i  =  1,  2,  3)  depends  on 
xt,  x2  only,  a  general  solution  to  (1.2)  is,  in  matrix  notation, 

u  =  a  f(z),  z  =  Xj  +  pxg.  (1.3) 

In  the  above  f  is  an  arbitrary  function  of  z,  and  p  and  a  are  determined 
by  inserting  (1.3)  into  (1.2).  We  have 
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{Q  +  p(R  +  RT)  +  p2T}  a  =  0  (1.4) 

where  the  superscript  T  denotes  the  transpose  and  Q,  R,  T  are  3x3  real 
matrices  whose  components  are 

^ik  =  Cilkl*  Rik  =  Cilk2*  Tik  =  Ci2k2‘  (LS) 

We  note  that  Q  and  T  are  symmetric  and,  subject  to  the  positiveness  of 
strain  energy,  positive  definite.  The  stresses  obtained  by  substituting 
(1.3)  into  (1.1)  can  be  written  in  terms  of  the  stress  function  ^  as 

Jil  =  ~^i,2’  ^i2  “  ^i,l»  (L6) 

in  which 

♦  =  b  £(*),  (1.7) 

b  =  (RT  +  pT)»  =  -  i(Q  +  pR)a.  (1.8) 

The  second  equality  in  (1.8)  follows  from  (1.4).  It  suffices  therefore  to 
consider  the  stress  function  $  because  the  stresses  ?ij  can  be  obtained  by 
differentiation. 

There  are  six  eigenvalues  p  and  six  eigenvectors  a  from  (1.4). 
Since  p  cannot  be  real  if  the  strain  energy  is  positive  [4],  there  axe  three 
pairs  of  complex  conjugates  for  p.  If  p^  atf  b^  (a  =  1,  2, ..,6)  are  the 

eigenvalues  and  the  associated  eigenvectors  we  let 

to  Pa  >  0.  P q-4-3  =  Ptf  »o+3  =  V  ba+3  =  t1'9) 

(a  =  1,  2,  3),  where  Im  stands  for  the  imaginary  part,  the  overbar 
denotes  the  complex  conjugate  and  ba  is  related  to  aft  through  (1.8). 

Assuming  that  the  pQ  are  distinct,  the  general  solutions  for  u  and  $ 

obtained  by  superposing  six  solutions  of  the  form  (1.3)  and  (1.7)  are 

3 

■  =  XKfc*(za>  +  »af  0+3^}' 
a=l 

3  (1-10) 

♦=XKw+Va«M- 

0=1 

In  (1.10)  fi,  fa,...,  fg  are  arbitrary  functions  of  their  argument  and 
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zcT  *1  +  Pa*2' 

The  above  fonnalism  is  due  to  Stroh  [1,2].  In  applications  all  we 
have  to  determine  is  the  form  of  the  arbitrary  functions  f  .  What 

distinguishes  the  Stroh  formalism  from  others  is  that  there  are  relations 
between  aQ  and  ba  which  allow  us  to  find  the  solution  easily  and/or  to 

simplify  the  solution  obtained.  These  relations  and  the  Barnett— Lothe 
integral  formalism  are  presented  next. 

In  closing  this  section  we  note  that,  in  most  applications,  fft  has 
the  same  function  form  sc  that  we  may  write 

W  =  4  </(*£.)' 

fo+3^o)  =  “  =  *'  2'  3’ 

where  are  arbitrary  constants.  The  second  equation  is  for  obtaining 
real  solutions  for  u  and  Equations  (1-10)  can  then  be  written  as 

3  3 

n  =  2  Re^V/'a*.  ♦  =  2  (111) 

a=l  a=l 


2.  THE  BARNETT— LOTHE  TENSORS.  The  two  equations  in  (1.8) 
can  be  rewritten  as 


where  I  is  the  3x3  identity  matrix.  Multiplying  both  sides  by  the  matrix 


leads  to  the  standard  eigenrelation  [5,6] 

=  p{,  (2.1) 


N1  N2 

N,  N* 


(2.2) 
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1 


—  x  j  —  m 


It  is  dear  that  Nj  and  Nj  are  symmetric  and  Nj  is  positive  definite.  It 
can  be  shown  that  — Nj  is  positive  semi— definite  [7].  Moreover,  — Nt  and 
— Nj  have  the  structure 


* 

* 

_ 4 

[*  0  *1 

N,= 

*  0  * 

.*  0  *. 

II 

r> 

T 

o  # 

o  o 

o  * 

■ 

(2.4) 


in  which  the  *  denotes  a  possibly  non-zero  element.  These  *  elements 
have  surprisingly  simple  expressions  in  terms  of  elastic  compliances  [7]. 
The  structure  of  N,,  Nj  shown  in  (2.4)  plays  important  roles  in  solving 
problems  and  interpreting  the  final  solutions. 

The  vector  (  =  (a,  b)  in  (2.2)  is  the  right  eigenvector  of  N.  It 
can  be  shown  that  (b,  a)  is  the  left  eigenvector.  The  left  and  right 
eigenvectors  associated  with  different  eigenvalues  are  orthogonal  to  each 
other.  The  orthogonality  relations  can  be  normalized  such  that 


Vb/J  +  ba'i/J={a0 


(2.5) 


where  Sn  is  the  Kronecker  delta.  Introducing  the  3x3  matrices  A  and  B 
by 


A  —  [aj,  a^,  a«j],  B  —  [b^,  bj],  (2.6) 


and  employing  (1.9),  the  orthogonality  relations  (2.5)  take  the  form 


bt  ati 

A  H 

I  0 

.  bt  at. 

B  B 

II 

.0  i. 

(2.7) 


The  two  6x6  matrices  on  the  left  are  the  inverses  of  each  other  and  their 
product  can  be  interchanged.  The  interchanged  product  is 

abt  +  XBt  =  i .  bat  +  3Xt, 

(2.8) 

aat  +  XXT  =  o  =  bbt  +  BBt. 

T 

Equations  (2.8)  tell  us  that  the  real  part  of  AB  is  1/2  and  that 
T  T 

AA  and  BB  are  purely  imaginary.  Hence  the  three  matrices  S,  H,  L, 
defined  by 

S  =  i(2ABT  —  I),  H  =  2iAAT,  L  = -2iBBT, 
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(2.9) 


are  real.  It  is  clear  that  H  and  L  are  symmetric.  It  can  be  shown  that 
they  are  positive  definite,  the  products 

SH,  LS,  ff^S,  SIT1 

are  antisymmetric,  and  the  relation 


holds  [3,6]. 


HL  -  SS  =  I 


(2.10) 


The  formulation  presented  so  far  assumes  that  the  eigenvalues  p 
are  distinct.  If  pi  =  P2,  say,  and  at  =  as,  the  solution  (1.10)  is  not 
general.  The  matrices  A  and  B  would  be  singular  and  the  orthogonality 
relation  (2.7)  is  not  valid.  Anisotropic  materials  for  which  pt  =  p2  and 
a(  =  as  are  called  degenerate  materials.  They  are  degenerate  in  the 
mathematical  sense,  not  necessarily  in  the  physical  sense.  Isotropic 
materials  are  a  special  case  of  degenerate  materials  for  which  pt  =  pa  = 
ps  =  i  and  a(  =  a?  *  a3.  In  many  applications  however  the  fink  solution 
depends  only  on  the  three  real  matrices  S,  H,  L  defined  in  (2.9).  Barnett 
and  Lothe  [3]  devised  an  integral  formalism  of  these  three  real  matrices 
which  circumvented  the  need  of  determining  the  eigenvalues  and  the 
eigenvectors.  Thus  the  problem  of  degenerate  materials  disappears.  The 
integral  formalism  is  as  follows.  Define  the  three  real  matrices 

Qik(fl>  =  cijk.Ys.  Rik(")  =  Tik(*>  =  cijk.mf,s,  (2  U> 

in  which  0  is  a  real  parameter  and 


ns  =  [cos0,  sin0,  0],  mj  =  [-sin0,  cos0,  0]. 

Equations  (2.11)  reduce  to  (1.5)  when  0=0.  Next  consider  the 
incomplete  integrals 


N1(»)  =  -T-1(*)RT(«),  n2(0)«  T-1(«), 

Nj(0)  =  R(0)T-1(0)RT(0)  -  Q(0). 
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where 


Ni(0)  reduce  to  Ni  in  (2.3)  when  0  =  0.  When  0  =  t  we  have  the 
complete  integrals  S(ir),  H(ir),  L(x).  Barnett  and  Lothe  proved  that  S, 
H,  L  of  (2.9)  are  identical  to  the  complete  integrals,  i.e., 

S  =  S(t),  H  =  H(x),  L=L(x).  (2.13) 

Thus  S,  H,  L  are  called  the  Bamett-Lothe  tensors  and  S(0),  H(0),  L(#l 
the  associated  tensors.  In  the  sequel,  dependence  of  S(0),  H(0),  L(0) 
on  0  will  be  given  explicitly  unless  0  =  ir,  and  dependence  of  Ni(0)  on  0 
will  be  given  explicitly  unless  0=0. 


As  we  see  from  the  integrals  in  (2.12),  there  is  no  need  to 
determine  the  eigenvalues  p  and  the  associated  eigenvectors  a  and  b. 
This  is  a  remarkable  result  which  has  been  widely  used  in  the  analysis  of 
anisotropic  elasticity.  It  should  be  pointed  out  that  there  are  cases  in 
which  the  final  solution  cannot  be  presented  entirely  in  terms  of 
Barnett— Lothe  tensors  and  their  associated  tensors.  In  that  case  we 
have  to  modify  the  general  solution  (1.10)  for  degenerate  materials  [8,9]. 


For  isotropic  elastic  materials  use  of  (2.12)  leads  to 


S  =  s 


0 

-1 

0 

.  H  =  - 

1 

0 

0 

0 

0 

0 

l-s» 

K 

0 

0 


0 

1— S3 
K 

0 


0 

0 

1 


K 

0 

0  ■ 

,  L  =  M 

0 

K 

0 

.0 

0 

1 

(2.14) 


where  /x  is  the  shear  modulus, 

1  .  _  1-2  v 

*  =  T=p  s  = 

and  v  is  the  Poisson  ratio.  For  general  anisotropic  materials  the 
structure  of  S,  H,  L  is  more  complicated.  For  orthotropic  materials  and 
for  monoclinic  materials  with  the  plane  of  symmetry  at  X3  =  0,  explicit 
expressions  of  S,  H,  L  are  obtained  in  [10,11].  We  will  show  later  that,  if 
a  proper  oblique  coordinate  system  is  chosen  as  the  natural  base  of  the 
tensors  S,  H,  L,  the  tensor  components  S'j,  HU  and  Lij  for  general 
anisotropic  materials  have  the  exact  expressions  as  that  shown  in  (2.14) 
for  isotropic  materials. 


3.  PHYSICAL  MEANINGS  OF  THE  EIGENVECTORS  a  AND  b.  Let 
a',  a"  be  the  real  and  imaginary  parts  of  a, 

a  =  a’  +  ia". 

A  complex  vector  is  also  called  a  bivector  [12,13].  The  real  vectors  a’ 
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and  a"  span  a  plane.  If  a  is  obtained  by  multiplying  a  by  a  complex 
factor  e*^  where  ip  is  real, 

a  =  e^a  =  a*  +  ia", 

in  which 


a’  =  a’cos^  -  a"sin  ip, 
a"  =  a’sin^  +  a"cos^- 


(3.1) 


Thus  the  real  and  imaginary  parts  of  a  lie  on  the  plane  spanned  by  a’ 
and  a".  Therefore  the  plane  is  called  the  polarization  plane  of  a,  or 
simply  the  plane  a,  which  is  invariant  with  the  multiplication  factor  on 

a.  As  ip  varies  (3.1)  show  that  a*  and  a"  trace  an  ellipse.  A  pair  of 
diameters  in  an  ellipse  is  said  to  be  conjugate  if  all  chords  parallel  to  one 
diameter  are  bisected  by  the  other  diameter.  Therefore  the  tangent  at 
the  extremity  of  one  diameter  is  parallel  to  the  other  diameter.  It  can  be 

shown  that  a’  and  a"  form  a  pair  of  conjugate  radii.  One  could  choose  a 

ip  such  that  a’  and  a"  are  orthogonal  and  hence  are  the  principal  radii  of 
the  ellipse  [14]. 

It  is  clear  that  the  bi vector  a  and  its  complex  conjugate  a 
define  the  same  polarization  plane. 

Consider  now  the  solution  (1.3).  The  displacement  u  is  a 
bivector  a  multiplied  by  f(z).  Regardless  of  the  position  (xi,  X2),  f(z)  is  a 

complex  factor  of  the  form  pe1^  where  p  is  real.  Whether  we  take  the 
real  or  imaginary  part  of  a  f(z),  u  is  polarized  on  the  plane  a  for  all 
(xi,  x2).  Likewise,  the  stress  function  #  of  (1.7)  is  polarized  on  the 
plane  b.  If  tp  is  the  surface  traction  vector  on  a  curved  boundary  T,  it 

can  be  shown  from  (1.6)  that 

*T  =  %  (3-2) 

where  rj  is  the  arclength  of  T  measured  in  the  direction  such  that  the 
material  is  located  on  the  right  hand  side  of  T.  Equations  (1.6)i  and 
(1.6)2  ue  special  cases  of  (3.2)  when  T  is  the  surface  xt  =  constant  and 
X2  =  constant,  respectively.  Since  ^  is  polarized  on  the  plane  b,  (3.2) 
tells  us  that  the  surface  traction  tp  is  polarized  on  the  plane  b. 

The  general  solution  (1.10)  or  (1.11)  implies  that  there  are  three 
polarization  planes  a(,  a2,  aj  for  the  displacement  u  and  three 
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polarization  planes  bt,  b2,  u3  for  the  surface  traction  tp.  For  monoclinic 

materials  with  the  plane  of  symmetry  at  x3  =  0,  at,  a2,  bt,  bj  all  define 
the  same  plane,  namely,  the  (xt,  x3)  plane.  As  to  a3  and  b3l  their  real 
and  imaginary  parts  are  parallel.  The  polarization  planes  degenerate 
into  lines  parallel  to  the  x 3— axis.  The  displacement  associated  with  a3 
and  the  surface  traction  tp  associated  with  b3  are  in  the  x3  direction. 

In  summary,  there  are  three  independent  (or  three 
one-component )  solutions  for  general  anisotropic  materials.  The 
displacement  ot  a  one-component  solution  is  polarized  on  the  plane  a 
while  the  surface  traction  on  any  boundary  is  polarized  on  the  plane  b. 
To  satisfy  a  prescribed  boundary  condition,  all  three  one-component 
solutions  are  in  general  needed.  In  surface  waves,  there  are 
one— component  surface  waves  (15,  16]  and  two— component  surface  waves 
[17,  18].  For  Green’s  functions  for  tne  infinite  space  due  to  a  line  force 
and  a  line  dislocation,  there  are  one-component  Green’s  functions.  The 
latter  will  be  discussed  in  Section  S. 


4.  THE  S  TENSOR.  Of  the  three  Barnett-Lothe  tensors,  the  tensor  S 
is  the  most  interesting  one.  By  writing  S  as 

S  =  ir'(LS),  (4.1) 

S  is  the  product  of  the  symmetric  positive  definite  tensor  h  ’  and  the 
antisymmetric  tensor  LS.  It  has  the  property  that 

tr  S  =  0,  det  S  =  0. 

Therefore  the  eigenvalues  of  S  are  0  and  ris  where 

s  =  {-^tr(S2)}1/2  (4.2) 

Denoting  the  associated  eigenvectors  by  e3  and  e^iej  where  ej,  e3,  e3 
are  real,  we  have 


S(et±ie3)  =  =Fis(ei±ie2),  Se3  =  0.  (4.3) 

Thus  e3  is  the  right  null  vector  of  S  and  ei±iej  are  the  right 

eigenvectors.  The  new  right  eigenvectors  e^iej  obtained  by 
multiplying  e^ie?  by  a  complex  factor  span  the  same  plane  as  e^iej. 
Therefore  the  plane  spanned  by  (et,  ej)  is  called  the  right  eigenplane. 

Let  e‘,  eJ,  e3  be  the  reciprocal  of  ei,  ej,  e3  so  that 

e‘.ej  =  #j. 

25 


(4.4) 


Consider  the  following  tensor  components  of  S,  H,  L  : 

S  =  S‘jei»ei,  H  =  Hiiei«ej,  L=Lijei«ei.  (4.5) 

Using  (2.10)  and  the  fact  that  SH,  LS  are  antisymmetric,  the  matrices 
formed  by  Sij,  HU,  Ljj  can  be  shown  to  have  the  structure  given  in  (2.14) 
where  s,  m,  *  are  now  independent  constants  [19).  Thus  as  far  as  the 
Barnett— Lothe  tensors  are  concerned,  anisotropic  materials  are  identical 
to  isotropic  materials  if  we  choose  an  oblique  coordinate  system 
represented  by  ei,  ej,  ej.  For  isotropic  materials  ei,  ej,  ej  are  unit 
vectors  in  the  direction  of  the  xj,  xj,  xj  axis,  respectively. 

It  should  be  pointed  out  that  (el,  e3)  and  e3  are,  respectively, 
the  left  eigenpiane  and  the  left  null  vector  of  S.  Do  ei,  ei  have  physical 
interpretations?  They  do.  They  are  explained  in  the  next  Section. 


5.  GREEN'S  FUNCTIONS  FOR  LINE  FORCES  AND  LINE 
DISLOCATIONS  IN  THE  INFINITE  SPACE.  There  are  several 
interesting  properties  associated  with  Green’s  functions  for  the  infinite 
space  due  to  a  line  force  /  and  a  line  dislocation  with  Burgers  vector  b 
applied  along  the  x3  axis.  The  basic  solution  is  obtained  from  (1.11)  by 
choosing  the  function  f(zQ)  such  that 

3  3 

u=4lmY  a  a  In  z  .  <fc  =  ^  Im  V  b  a  In  z  .  (5.1) 

r  ^  cfa  a'  T  r  ^  a*a  a  v  ’ 

a=l  a=l 


Since  In  zft  is  a  multi-valued  function  we  introduce  a  cut  along  the 
negative  xp- axis.  In  the  polar  coordinate  system 

Xj  =  r  cosd,  x^  =  r  sin0,  (5.2) 

the  solution  (5.1)  applies  to 


Therefore 


— t  <  0  <  x,  r  >  0. 


In  za  =  In  r  ±  iir  at  0  =  ±t,  for  a  =  1,  2,  3.  (5.3) 

Equations  (5.1)  represent  three  one-component  Green’s 
functions.  For  each  a,  u  is  polarized  on  the  plane  aQ  and  the  surface 

traction  tp  is  polarized  on  the  plane  b^.  The  discontinuities  in  u  and  $ 

across  $  =  *t  are,  respectively,  the  line  dislocation  ba  and  the  line  force 
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fa  for  the  one— component  Green’s  function.  Hence  by  (5.3), 

ba  =  2  Ret^qJ,  /“  =  2  Re(baqo),  (5.4) 

which  show  that  ba  is  on  the  plane  aft  and  /  a  is  on  the  plane  b^.  We 
therefore  have  the  result  that  the  one-component  Green’s  function  has  n 
and  ba  polarized  on  the  plane  a  and  has  fa  and  the  surface  traction  tp 
polarized  on  the  plane  b^. 

To  obtain  a  one-component  Green’s  function  we  may  assume  an 
arbitrary  complex  constant  q^.  Equations  (5.4)  then  provide  ba  and  fa 
required  for  the  one— component  Green’s  function.  Alternately  we  may 
prescribe  an  /  a  which  lies  on  the  plane  bft.  Equation  (5.4)2  can  be 

solved  for  qfl  and  (5.4),  gives  the  associated  ba.  To  solve  (5.4)2  for  qtf 
let  the  real  and  imaginary  parts  of  ba  and  qa  be  written  as 


We  then  have 


b  =  b’  +  ib”  .  q  =  q’  +  iq" . 
a  a  a  ’  Ha  Ha  4o 


/“=2(>>X-W 


from  which  q  ’  and  q"  can  be  determined. 

’a  a 

When  /  and  b  are  prescribed  arbitrarily,  we  need  all  three 
one-component  Green’s  functions  for  the  solution.  Making  use  of  (2.6), 
(5.1)  are  rewritten  as 

u  =  ^  Im{  A<ln  z>q},  $  =  ~  Im{B<ln  z>q),  (5.5) 

in  which 

T 

q  =  [qi.  qa,  qa] 

and 

<ln  z>  =  diag[ln  zj,  In  Z2,  In  Z3] 


is  a  diagonal  matrix.  Equations  (5.5)  must  satisfy  the  conditions 


u(x)  -  u(— ir)  =  b, 


«*)-«-*)  =  f, 

27 


which  lead  to 


2  Re(Aq)  =  b,  2  Re(Bq)  =  / 
This  can  be  written  as 


It  follows  from  (2.7)  that 


bt  at 
bt  at 


Hence 


or 


q  —  AT/+  BT6, 
1a  =  V-f+  bal 


(5.7a) 

(5.7b) 


Inserting  (5.7b)  into  (5.4)  gives  us  ba  and  fa  in  terms  of  b  and  / 

We  show  next  that  the  solution  (5.5)  together  with  f5.7a)  can  be 
rewritten  into  a  real  form.  Equations  (2.9)  are  identities  wnich  convert 
certain  combinations  of  complex  quantities  involving  A,  B  to  real 
quantities  S,  H  and  L.  The  following  identities  are  useful  for  problems 
related  to  line  forces  and  line  dislocations  [20]. 

2  Acln  z>AT  =  — i[(ln  r)I  +  *S(0)]H  +  *H(0)[I  —  iST], 

2  B<ln  z>  ’*•  =  [(In  r)I  +  *ST(0)](I  -  iST]  +  iirL(0), 

(5.8) 

2  Acln  z>BT  =  [(In  r)I  +  ttS(^)][I  —  iS]  +  i*H(0)L, 

2  Bcln  z>BT  =  i[(ln  r)I  +  tST(^)]L  -  irL(0)[I  -  iS). 


These  identities  allow  us  to  convert  the  complex  expressions  on  the  left 
to  real  quantities  shown  on  the  right  which  are  obtainable  directly  in 
terms  of  elasticity  constants  through  (2.12)  and  (2.13).  With  the 
identities  (5.8),  the  solution  (5.5)  together  with  (5.7a)  is  converted  into  a 
real  form  as 
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(5.9) 


2u  =  -  4  (In  r)h  -  S(9)h  +  H(9)g, 

7T 

2+  =  i(lnr)g+L(9)h+ST(«)g, 

where 

g  =  Li-  ST/  h  =  S6  +  H/  (5.10) 

From  (3.2),  the  surface  traction  tg  on  any  radial  plane  0  =  constant  is  in 

the  direction  of  g  which  is  invariant  with  the  choice  of  the  radial  plane. 
The  infinite  displacement  no  at  r  ==  0  is  in  the  direction  of  h.  Moreover, 
the  relation  [14] 

g-h=  f-b 

is  easily  established  using  (2.10)  and  the  anti-symmetric  property  of  LS 
and  SH. 


We  now  present  physical  interpretations  of  ei  and  eh  Using 
(4.5)  and  the  discussions  following  (5.10),  (5.10)  can  be  written  as 


2*rt0  =  g  =  [Lij(ei  ♦  6)  -  Sii(ej-/)]e*, 

-2r(ln  r)-1u0  =  h  =  [S^ei*  6)  +  Hii(eJ*/)]ei. 

With  the  structure  of  S»j,  HU,  Ly  shown  in  (2.14)  and  using  (4.4),  it  can 
be  shown  that  if  6  is  along  e3  and  /  is  along  e3,  no  is  in  the  direction  of  e$ 
and  t^  in  the  direction  of  e3.  On  the  other  hand,  if  b  is  on  the  right 

eigenplane  (et,  e2)  and  /  is  on  the  left  eigenplane  (el,  e3),  no  is  polarized 
on  the  right  eigenplane  and  t^  is  polarized  on  the  left  eigenplane.  More 

relations  between  ei,  e2,  el,eJ  and  b  and  /  can  be  found  in  [14]. 


6.  BIMATERIALS  AND  INTERFACE  CRACKS.  Let  0  =  0O  be  the 
interface  between  two  materials  in  the  bimaterial.  The  half-space  9o  < 
6  <  0o+*  is  occupied  by  material  1  and  the  other  half-space  0<r- r  <  6  < 
0o  is  occupied  by  material  2.  They  are  rigidly  bonded  together  along  0  = 
0o-  For  a  line  force  /  and  a  line  dislocation  I  applied  at  the  origin  r  =  0, 
(5.9)  is  a  basic  solution  which  applies  to  both  materials.  We  may  add 
constant  terms  to  the  right  hand  sides  of  (5.9)  which  produce  a  rigid 
body  displacement  but  no  stresses.  Therefore  consider  the  solution 

2  o,(r, 9)  =  -  i(In  r)h  -  [S,(9)-S,(90)lh  +  (H,( 

2  ♦.(!,*)  =  j  (In  t)g  +  (L,(9)— L,(9o))h  +(sT(«)-sT(«o)lg, 
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(6.1) 


for  material  1  in  do  <  0  <  9q+t.  The  subscript  1  denotes  material  1. 
The  solution  for  material  2  is  obtained  from  (6.1)  by  replacing  the 
subscript  1  by  2  while  keeping  the  same  constants  g  and  h.  It  is  readily 
shown  that  the  continuity  of  n  and  $  at  0  =  0O  is  automatically  satisfied. 
The  discontinuity  in  u  and  $  across  d  =  d0  ±  r,  which  should  be  equal  to 
6  and  / ,  respectively,  leads  to  two  equations  for  g  and  h  which  are 
independent  of  d0  [21].  Therefore,  the  stresses  obtained  by  substituting 
$i  of  (6.1)  and  similar  equation  for  ^2  into  (1.61  are  independent  of  the 
location  d0  of  the  interface  !  This  unexpert  ea  phenomenon  defies  an 
intuitive  explanation  even  for  isotropic  bimaterials. 

One  of  the  most  studied  problems  in  anisotropic  elasticity  is  the 
problem  of  interface  cracks  in  bimaterials  [22-33].  Let  X2  >  0  be 
occupied  by  material  1  and  x2  <  0  be  occupied  by  material  2.  The 
interface  crack  of  length  2  a  is  located  at 

x2  =  0,  |x,|  <  a. 

The  bimaterial  is  subject  to  a  uniform  traction  tp  and  — tp  at  the  crack 

surfaces  x2  =  +0  and  -0,  respectively.  The  stress  singularities  near  a 

tip  of  the  interface  crack  is  proportional  to  r^  where  r  is  the  radial 
distance  from  the  crack  tip  and  6  is  a  constant  depending  on  the  material 
property  of  the  bimaterial.  It  is  shown  in  [24]  that  there  are  three 
singularities  given  by 

<*=  -5,  “j  +  i%  and  -\-H, 

where 

7  =  EF  lnlH? =  x tanh  1/3’ 

0=  [-jtr(S2)]1/2<l. 

In  the  above 

S  =  D-1W, 

D  =  L-‘  +  t  \  W  =  SlL-‘  -  Sjl^1, 

in  which  D  is  symmetric,  positive  definite  and  W  is  anti-symmetric. 

Thus  S  has  the  same  properties  as  the  S  tensor.  The  eigenvalues  of  S  are 
vi/9  and  0  and  the  associated  right  eigenvectors  are  denoted  by  di±id2  and 
do,  respectively.  The  left  eigenvectors  can  be  shown  to  be  D(dividj)  and 
Ddo  [34].  Hence  d0l  Ddo  are,  respectively,  the  right  and  left  null  vectors 
while  the  planes  spanned  by  (di,  da)  and  (Ddi,  Dda)  are  the  right  and  left 
eigenplanes. 


(6.2) 

(6-3) 
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The  two  materials  in  the  bimaterial  are  said  to  be  mismatched 
when  W  *  0.  W  =  0  if  and  only  if  0  =  0  (and  hence  7=0)  [24,  251. 
For  mismatched  bimaterials  (7  f  0),  the  displacement  at  the  crack 
surface  is  oscillatory.  This  leads  to  the  physically  unacceptable 
interpenetration  of  the  crack  surfaces. 


When  0  =  0  the  solution  in  materials  1  and  2  both  have  the 
expression 


in  which 


u  =  Re  {A<f(z)>B_1}tr, 
4  =  Re  {B<f(z)>B—1}tp, 


(6.4) 


f(z) = y  za-o2 


z. 


Of  course  A,  B  and  z  in  material  1  and  material  2  would  be  different. 
There  is  no  oscillation  in  displacement  and  the  stress  is  the  square  root 
singularities. 


The  following  results  are  taken  from  [34].  When  0  t  0,  the 
solution  is  still  given  by  (6.4)  if  the  applied  traction  tp  is  the  null  vector 

of  W,  i.e.,  if 


or,  by  (6.3), 


Wtp  =  0, 

Stp  =  0. 


Thus  when  the  applied  traction  is  in  the  direction  of  the  right  null  vector 
do,  there  is  no  oscillation  in  displacement.  The  crack  surface  opening 

an  =  u(x t,  +0)  -  n(x,,  -0),  |x,|  <  a, 

is  in  the  direction  of  the  left  null  vector  Ddo  and  the  surface  traction  on 
the  surface  x2  =  0  outside  the  crack  is  in  the  direction  of  the  right  null 
vector  do. 


If  the  applied  traction  is  not  in  the  direction  of  do,  we 
decompose  it  into  two  components.  One  is  along  the  right  null  vector  do 
and  the  other  is  on  the  right  eigenplane  (d|,  dj).  Explicit  solutions 
associated  with  the  one  on  the  right  eigenplane  can  be  found  in  [34]  in 
which  the  displacement  is  oscillatory.  It  suffices  to  mention  that  the 

crack  surface  opening  au  lies  on  the  left  eigenplane  of  S  while  the  surface 

traction  along  the  surface  x2  =  0  lies  on  the  right  eigenplane  of  S. 
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DISCUSSION 


We  have  shown  that,  in  many  respects,  anisotropic  elastic 
materials  have  properties  which  are  similar  to,  or  generalization  of,  the 
properties  of  isotropic  materials.  Analogous  to  the  antiplane 
deformations  of  isotropic  materials,  anisotropic  materials  have 
deformations  which  are  polarized  in  one  direction  while  the  surface 
traction  vector  on  any  boundary  is  polarized  on  a  different  direction. 
Similar  to  the  in-plane  deformations  of  isotropic  materials,  anisotropic 
materials  can  have  deformations  which  are  polarized  on  one  oblir 
plane  while  the  surface  traction  vector  on  any  boundary  is  polarizer 
another  oblique  plane. 

Simple  problems  for  which  we  thought  we  have  understood  them 
thoroughly  still  yield  new  information  due  to  the  simplification  of  the 
solutions  by  the  Stroh  formalism.  For  example,  consider  the  Griffith 
crack  of  length  2 a  located  at  X2  =  0,  |xt|  <  a  in  the  infinite  anisotropic 
elastic  medium.  When  the  traction  applied  at  the  crack  surfaces  is  in 
the  direction  of  the  xj  axis,  the  crack  opening  is  in  general  not 
symmetric  with  the  x3  axis  as  expected.  However,  the  xi  axis  outside  the 
crack  remains  a  straight  line  (i.  e.,  the  uj  component  of  the  displacement 
along  the  xt  axis  vanishes).  If  the  traction  applied  at  the  crack  surface  is 

the  null  vector  of  SLT1,  all  three  displacement  components  along  the  xt 
axis  vanish.  If  the  applied  traction  is  in  the  direction  of  the  vector 
formed  from  the  second  column  of  L,  the  hoop  stress  vector  along  the 
crack  surface  is  independent  of  xy  (21j. 

Other  interesting  properties  worth  mentioning  are  the  physical 
implications  of  the  eigenvalues  p.  For  the  Green’s  functions  for  a 
half-space  subject  to  a  singularity  in  the  form  of  line  forces  and  line 
dislocations,  the  solution  can  be  obtained  by  a  superposition  of  the 
Green’s  function  due  to  the  same  singularity  for  the  infinite  space  and 
several  image  singularities  located  outside  of  the  space  occupied  by  the 
material.  The  locations  of  the  image  singularities  are  determined 
exclusively  by  the  eigenvalues  p.  Moreover,  the  locations  of  the  image 
singularities  are  independent  of  the  nature  of  the  singularities  concerned 
[35j.  If  the  singularities  are  line  forces  and  line  dislocations,  the  image 
singularities  are  also  line  forces  and  line  dislocations.  For  degenerate 
materials  for  which  isotropic  materials  are  a  special  case,  two  or  more  of 
the  image  singularities  coalesce  into  one  singularity,  creating  a  new 
singularity  in  the  form  of  a  double  force,  a  concentrated  couple,  and/or  a 
higher  order  singularity  which  are  well  known  for  isotropic  materials  [36] 
but  have  not  been  satisfactorily  explained  in  the  past. 
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